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Abstract 



> 

' We critically review some concepts underlying current applications of gravity 

Q . theories with Lagrangians L = f{g^u,Ra/3fiu) to cosmology to account for the 

CN I accelerated expansion of the universe. We argue that one cannot reconstruct the 

function / from astronomical observations either in cosmology or in the solar 
. system. The Robertson- Walker spacetime is so simple and "flexible" that any 

Q I cosmic evolution may be htted by infinite number of various Lagrangians. We 

O^' show on the example of Newton's gravity that one cannot recover the correct 

5-H . equation of motion from its approximate solution. Any gravity theory different 

from general relativity generates a new cosmological theory and all the successes of 
the standard cosmological model are lost even if a single solution of the theory well 
^ ■ fits the observations. Prior to application of a given gravity theory to cosmology or 

elsewhere it is necessary to establish its physical contents and viability. This study 
may be performed by a universal method of Legendre transforming the initial 
Lagrangian in a Helmholtz Lagrangian. In this formalism Lagrange equations of 
motion are of second order and are the Einstein field equations with additional 
massive spin-zero and spin-two fields. All the gravity theories differ only by a 
form of interaction terms of the two fields and the metric. Initial conditions for 
the two fields in the gravitational triplet depend on which frame (i.e., the set 
of dynamical variables) is physical (i.e. matter is minimally coupled in it). This 
fact and the multiplicity of possible frames obstruct confrontation of solutions 
to equations of motion with the observational data. A fundamental and easily 
applicable criterion of viability of any gravity theory is the existence of a stable 
ground state solution being either Minkowski, de Sitter or anti-de Sitter space. 
Stability of the ground state is independent of which frame is physical. 
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PACS numbers: 04.50. +h, 98.80.Jk 



1 Introduction: problems with gravity theories 

Metric theories of gravity, where the Lagrangian is any smooth scalar function of the 
curvature tensor, L = f{g^u,Rai3nv)^ named metric nonhnear gravity (NLG) theories, 
have first attracted attention as possible candidates for foundations of quantum gravity 
due to their renormalizability properties [Ij and then as a possible source of inflationary 
evolution of the very early universe [2] . A recent revival of interest in these theories has 
come from cosmology. In fact, the theoretical state of affairs in cosmology is astonishing. 
The universe consists of baryons (4 percent), unknown stable massive particles forming 
nonbaryonic dark matter which do not fit the standard model of particle physics (26 
percent) and dark energy (70 percent) about which we have only negative knowledge: 
these are not particles. About 96 percent of the material content of the universe is a 
great mysterious puzzle. 

On the other hand the laboratory experiments and astronomical observations con- 
firming general relativity are still not very numerous and belong to a rather narrow class 
of tests. It is therefore attractive to conjecture that both the gravitational stability of 
galaxy clusters and the acceleration of the universe are not due to some unknown forms 
of matter, but can be accounted for by some modification of gravity theory. Modifi- 
cations may go in all possible directions; here we deal with the most popular concept, 
the metric NLG theories. Among these the restricted NLG theories, wherein the La- 
grangian is a function of the curvature scalar alone, L = f{R), have been most frequently 
investigated. This approach to the dark matter and dark energy problem is sometimes 
referred to as a "curvature quintessence scenario". A typical motivation underlying this 
approach is following. Consider a Lagrangian of the form L = R + R"^ + ^ and the 
Robertson-Walker (R-W) spacetime. In the very early universe, when the curvature 
was large, the R^ term was dominating generating some kind of inflation |2j . At present 
R is small and the term dominates giving hopefully rise to the accelerated evolution. 
And for the most time in the history the curvature scalar had intermediate values so 
that the linear term was leading preserving all the successes of the standard Friedmann 
cosmology [SI IH El El El El E] • This argument, attractive as it sounds, is misleading for 
three reasons. 

Firstly, if one makes a correction to the Einstein-Hilbert Lagrangian in the form 
L = R + £{R) where £{R) is any nonlinear function, it is not true that the resulting 
corrections to solutions of Einstein's theory are small when e{R) is very small and 
become significant only when e{R) is sufficiently large. The point is that any nonlinear 
correction to L = R drastically alters the dynamical structure of the theory: the field 
equations become of fourth-order instead of second order and the higher curvature 
terms, even seemingly small, are always very important. To show the effect we consider 
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a very simple model, a one-dimensional harmonic oscillator perturbed by introducing a 
small term with the third derivative, 



e X +x + LO^x = 0, 



with l^l -C 1. One seeks for solutions of the form x — e , then A = const is a solution 
of a cubic equation eX^ + + o;^ = 0. One solves it by perturbing the two unperturbed 
solutions, i.e. one sets X± = ±iuj + ea±. Up to terms linear in e the solutions of the 



being slowly damped or amplified periodic solutions for the unperturbed oscillator. 
These are approximations to exact solutions which are analytic in e at £ = 0. However 
there is also a non-analytic at e = solution of the differential equation corresponding 
to a third root of the cubic equation. Assuming that A3 is of order and keeping only 
the leading terms (of order £~^) one gets A3 = — ^ and the third solution is exponentially 
growing {e < 0) or fading. 



This solution is qualitatively different from the other two and exists for arbitrarily small 
e. One expects that the space of solutions is significantly extended by qualitatively new 
solutions due to any correction e{R) to L — R. 

Secondly, the Lagrangian is not a physical observable whose value or variability 
gives some insight into the state of a given physical system. It is a kind of a generating 
function giving rise to equations of motion for the system and observables such as 
energy-momentum tensors (canonical and variational). Any gravitational Lagrangian 
is, by definition, made up of scalars, in general these are all invariants of the Riemann 
tensor (and possibly their derivatives), while the resulting field equations are tensor 
ones. Any assumption about the value of a scalar appearing in L, say R, actually 
tells very little about corresponding solutions. For instance, setting R = in general 
relativity one gets not only all vacuum solutions but also those for matter with a traceless 
energy-momentum tensor. Actually one can say something nontrivial about the sought 
for solution merely by inspection of the Lagrangian L = f{R) only in the case of the 
simplest non-maximally symmetric spacetime, the Robertson- Walker one. In fact, the 
Riemann tensor for this metric is determined by the cosmic scale factor a{t) and for 
any gravitational Lagrangian (also that explicitly depending on the Weyl curvature) the 
field equations reduce to one quasilinear third order ODE for a{t), hereafter named the 
quasi-Friedmannian equation. Then assuming that it! is large in some epoch of cosmic 
evolution and small in another one, one may neglect small terms in this equation and 
find approximate (or even exact) solutions in these epochs. In so doing one must take 



cubic equation are a± = 



Thus one has two almost periodic solutions 




X3 = a3exp( — ). 
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care since there is no easily available information about values of R in various eras of the 
cosmic evolution. In fact, in general relativity R = —StcGT where T is the trace of the 
matter energy-momentum tensor and for cosmic perfect fluid T = 3p — p. In the early 
universe the ultrarelativistic plasma has the equation of state p = ■^p and then R = 0. 
Hence in the transition from the radiation to the matter era R increases, contrary to 
what has usually been assumed in the motivation for L = R + R^ + ^ and only later on 
it slowly fades. What is more important, the relationship between R and T (as well as 
between and T^j^) is algebraic only in general relativity and for a nonlinear L = f{R) 
it becomes a differential equation for R. If the field equations (in Jordan frame) are 
Efj,u{g) = SttGT^u, then their trace reads 

Rf'{R) - 2f{R) + 3f"'{R)g^"'R,f,R,^ + 3f"{R)DR = StxGT 

and one cannot a priori estimate which terms in the equation are dominant and which 
ones may be neglected. To this end one must assume not only an approximate value of 
i? in a given epoch but also that R is slowly and monotonically varying in order to esti- 
mate the derivatives; in this way one eliminates possible rapidly oscillating solutions. All 
this makes sense if one is convinced that he deals with the correct Lagrangian, e.g. the 
Lagrangian has been derived from first principles (string theory, quantum gravity etc.) 
or otherwise motivated. Going in the opposite direction, i.e. attempting, as is recently 
done, to reconstruct the underlying Lagrangian from observed qualitative features of 
the cosmic scale factor means that one has to construct the whole relativistic cosmology 
anew. 

In fact, from the general Hawking-Penrose singularity theorem, valid in general rel- 
ativity, it follows that our universe contains a singularity since the cosmic fluid satisfies 
the strong energy condition. Then in the case of the Robert son- Walker spacetime the 
Friedmann equation implies that the singularity was in the past, the initial Big Bang, 
and the cosmic scale factor monotonically grows from zero at the curvature singularity. 
We stress that these are generic quahtative properties of any solution to the Friedmann 
equation, independent of a specific equation of state for the cosmic fluid matter. On the 
other hand for a generic L = f{R) the singularity theorem cannot hold. Whether or not 
the theorem holds must be proved case by casqj. For an arbitrary function f{R), even 
for many of those Lagrangians which admit a solution qualitatively fitting the astro- 
nomical data (an acceleration phase at present preceded by a deceleration phase), there 
was no initial singularity. There was no Big Bang and the following "early universe" 
when it was small, dense and hot. And the quasi-Friedmannian equation for such a 
Lagrangian does not ensure that the cosmic scale factor grows monotonically from an 
initial small value. What is actually done by those authors who say "today the R~^ term 

^For theories in which the theorem holds it takes a form different from that in general relativity 
since R^^ is not a function of T^^. 
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is leading in L while in the early universe the term was dominant" is merely choos- 
ing a class of (approximate) solutions for which a and R significantly and monotonically 
vary in the course of a cosmic evolution; the solutions in this class are, by construction, 
ever increasing. This is a kind of "fine tuning" since one may expect that there are 
classes of qualitatively different solutions which are no less typical. This conjecture is 
at least partially supported by rigorous investigations in the Einstein frame where the 
third gravitational degree of freedom for L = f{R) gravity is revealed as a nonlinear 
scalar field minimally coupled to Einstein gravity: for certain simple (power-law) scalar 
field potentials all solutions (in the spatially fiat Robertson- Walker spacetime) are os- 
cillatory at late times while for potentials which are bounded from above there exist 
solutions which are global to the past with the Hubble parameter H = a/a converging 
to a constant nonzero value as t — > —oo (TU]. Taking into account the complexity of the 
quasi-Friedmannian equation for any nonlinear f{R), investigation of the qualitative 
properties of all solutions is not easy and can only be done (case by case) by performing 
a phase-space analysis. Is is fortunate that the equation can be reduced to a first order 
equation for H viewed as a function of H in the case of the Einstein-de Sitter spacetime 
(k = 0) and a second order one for the open and closed univers^ (Starobinsky in ref. 
[2]). A preliminary analysis of evolution of the Einstein-de Sitter metric in the case od 
a couple of special Lagrangians was performed by Carroll et al. [1]. For the Lagrangian 
most frequently studied, introduced in [3], they found an attractor solution a — > for 
late times and an exact solution a oc t^/^ starting from a curvature singularity at t = 0. 
However it is difficult to see from their phase portraits, which cover only a piece of one 
quadrant of the phase space, whether all solutions emerge from the singularity or there 
are nonsingular solutions (besides the exponentially growing and decreasing ones with 
constant H) and whether oscillatory (non-monotonic) solutions are excluded. It should 
be emphasized that if for a given accepted Lagrangian the two questions are answered 
"no", then all the successes of the standard Friedmann cosmology are lost. 

The last statement leads us to the third reason, taking a form of a problem: to what 
extent can one reconstruct the Lagrangian from a given solution? Clearly any given 
function may be viewed as a solution of many diverse differential equations. Requirement 
that the function is a solution of a Lagrange equation imposes stringent restrictions on 
possible equations and allows one to look for a unique answer. A simple example shows 
that there are cases when it can be effectively done under some conditions. The Newton's 
gravitational force may be expressed as the gradient of a potential and one seeks for a 
corresponding Lagrangian. The potential generates a differential scalar S = di(f) dicj) {xi 
are the Cartesian coordinates) and the kinetic part of any Lagrangian should be some 
function f{S). Furthermore there may be a potential energy V^(0) and a candidate 

^There is a subtle mathematical assumption necessary to decrease the order of the equation, namely 
that Hit) is a monotonic function. Then the analysis loses all solutions ait) which are not monotonic. 
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Lagrangian is L = f{S) — V{(j)). Inserting the Newton's potential = — - into the 
Lagrange equation 

4/" + 2/'(S) A0 + V'i(t>) = 

= one gets 

-^nS) + V'{<P) = 0. 

The problem is indeterminate since there are two unknown functions and one equation. 
A unique solution corresponding to the Laplace equation arises either upon setting = 
or requiring linearity of the resulting equation. 

Suppose now one is attempting to determine the Lagrangian of an underlying restricted 
NLG theory from observations in the solar system. It is reasonable to assume that the 
metric is static spherically symmetric. As we discuss further in this work one of essential 
requirements that any viable L = f{R) theory should satisfy is that it admits solutions 
either with R^^i, = or i?^,^ = Ag^^, with A constant positive or negative (Einstein 
spaces); clearly there are infinitely many Lagrangians having this feature. Thus each 
L = f{R) theory under consideration has either Schwarzschild or Schwarzschild-(anti)de 
Sitter spacetime as a solution and if it is this solution that is realized in the nature then 
no set of observations and experiments can reconstruct the function f{R)- In these theo- 
ries the Birkhoff theorem does not hold and other static spherically symmetric solutions 
do exist (very few of them are known [TT]). Any solution different from R^^, = Kg^y 
is generated by nonzero values of a scalar field being a nongeometric component of a 
gravitational doublet of fields, see section 3. There are two possibilities depending on 
whether the doublet is in the Jordan frame (more precisely, in Helmholtz- Jordan frame) 
or in the Einstein frame. If the measurable quantities form the Jordan frame, then the 
spin-0 gravity has ordinary matter as a source. In the solar system the source (the sun) 
enforces spherical symmetry of the scalar. The field is massive and unless its mass is 
extremely small it is a short range one. The presence of the scalar gives rise to a non- 
Schwarzschildean solution which is very close to ordinary Schwarzschild one (possibly 
with A 7^ 0). Nevertheless a recently found approximate solution for a specific form of 
f{R) shows that even small corrections are detectabl^ and in this case they are ruled 
out by measurements [12]. In the case the Einstein frame is observable the scalar is 
independent of any local matter distribution and in particular in the solar system it is 
not determined by position and mass of the sun. It is rather of cosmological origin, 
i.e. the initial data are fixed and common for the entire universe. Then the scalar field 
is (almost) homogeneous and to avoid a conflict with the cosmological observations it 
must be very weak giving rise to unmeasurable effects in the solar system. This means 



•^This solution belongs to the class of solutions which are generated by the higher-than-second order 
of the field equations in JF, i.e. it is analogous to the third solution for the perturbed harmonic oscillator 
discussed above. The solution does not reduce to a Schwarzschildean one in the limit f{R) R and 
this is why it considerably deviates from the former even if the difference f{R) — R is negligibly small. 
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that all observations performed there will confirm the Schwarzschild metric. Even the 
cosmological constant is undetectable locally since the most sensitive to it effect, the 
perihelion shift of Mercury, requires |A| > 10~^^m~^ [13], while the cosmologically ad- 
missible value is |A| ^ 10~^^ m~^. Thus observations may effectively distinguish between 
the two frames. If Jordan frame is physical then small corrections to Schwarzschild met- 
ric generated by the scalar are ruled out by the observations in vicinity of the sun in the 
case of most Lagrangians investigated up to now [TU [TJ [T^l [HI [El [HI [IHl [IH], though 
there are arguments that the linear approximation applied to calculate these corrections 
may not be valid within the solar system |20j . 

The Robert son- Walker spacetime has a higher (six-dimensional) symmetry group 
and its metric contains an arbitrary function which may be freely adjusted so that each 
metric NLG theory, i.e. any L = f{g^u,Rai3^iv)^ admits it as a solution of the corre- 
sponding quasi-Friedmannian equation and in this sense this spacetime is "flexible" 
and universal while Minkowski, de Sitter and anti-de Sitter spaces are not. Hence one 
expects that infinitely many of the theories predict a deceleration phase in the past 
and an accelerated evolution in the present epoch. It has been conjectured that any 
cosmological evolution may be realized by some specific L = f{R) [H]- And it is im- 
pressive that the astronomical data may be fitted by such diverse functions as rational 
and exponential ones [21], a combination of two confluent hypergeometric functions [6], 
a combination of two hypergeometric functions [22j and even an implicit form of f{R) 
was found which is consistent with the three years collection of WMAP data [6j . Also a 
"toy model" based on a minimal curvature conjecture {R is always above some minimal 
value Rq > 0) fits exceptionally well the ACDM cosmological model [23]. These results 
were found by studying the quasi-Friedmannian equation, i.e. the dynamics of L = f{R) 
theory; if one applies the dynamically equivalent Einstein frame it is easy to show that 
for any scale factor a{t) there exists at least one function f{R) having this factor as a 
solution of the field equatioE0 Plj . Therefore there is no surprise that also by inclusion 
of corrections in a form of inverse squares of the Ricci and Riemann tensors the Su- 
pernovae data can be fitted without the need of any dark energy [261 H]. On the other 
hand in the metric-affine formalism of gravity theories it was found that the current 
cosmological data (the supernovae la "gold set", the CMBR shift parameter and the 
linear growth factor) impose stringent restrictions on the test Lagrangian L = R + aR^ 
since /3 ~ 0,05 [27J; this result, however, does not uniquely confirm the ACDM model 
and is rather due to the specific choice of the Lagrangian. 

The main problem of cosmology which arises in this field of investigations is therefore 

''Yet in a recent paper it is found that for a large class of the Lagrangians the evolution of linear 
density fluctuations does not fit the structure formation on both large and small scales simultaneously 
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whether it is possible at all to effectively uniquely recover the underlying Lagrangian 
L = f{R) from an exact analytic solution a{t) to the quasi-Friedmannian equation. 
(Reconstruction of a Lagrangian depending on more than one curvature invariant is 
impossible.) The problem was investigated in more detailed way in two papers. First 
Capozziello et al. in an ingenious work [28] attempted to determine f{R) assuming that 
one knows from observations a precise analytic dependence of the Hubble parameter H 
on the redshift parameter z, H{z). They expressed the curvature scalar i? as a func- 
tion of H{z) and then derived an equation for f[R{z)] from the quasi-Friedmannian 
equatioE@. The equation for / is a linear third order ODE with coefficients being ex- 
tremely complicated functions of H{z) and hence, in spite of its linearity, is analytically 
intractable. The authors were able to find out the initial data from observations for the 
equation so that there is a unique solution, but it is inaccessible. The method is unique 
and ineffective. It requires an exact analytic form of H{z) while astronomic observations 
provide only a finite set of values H{z) which are affected by large errors. To circum- 
vent this difficulty the authors used a couple of simple functions approximating the real 
dependence of H on z which arise in recent models for dark energy such as quintessence 
and the Chaplygin gas. Given these analytic expressions, the numerical integration of 
the equation for f[R{z)] is inevitable, then the solution is approximated by an empirical 
function. The final analytic fit for f{R) looks rather unconvincing (eq. (48) of [28] ) 
and the authors caution the reader against drawing any physical implications from it. 
Recently Fay et al. |29j have attempted to search for the form of f{R) corresponding 
to the ACDM cosmology applying the general autonomous system for the k = R-W 
spacetime. They found, among other functions, L = R — 2A, but this Lagrangian cor- 
responds to saddle points of the autonomous system while the attractors of the system 
do not reproduce the radiation and the matter eras of the universe and give rise to a 
complicated irrational function f{R)- Their work clearly shows that the spatially fiat 
R-W spacetime is particularly deceptive in this respect. In fact, the dynamics of this 
spacetime is not so simple as is usually assumed. We almost never employ the exact 
solutions of the Friedmann equation for the real matter content of the universe as they 
are too involved and we are so accustomed to the approximate solutions (which very well 
fit the observations) that we view them as exact ones. For example, the period when the 
energy densities of the "radiation" (all the ultrarelativistic particles including neutrinos) 
and of the "dust matter" (nonrelativistic baryons and possibly the dark matter) were 
comparable was longer than the whole radiation era and in this period a{t) was very 
complicated, nevertheless it is commonly assumed that in the radiation era there was no 
dust and this era was suddenly transformed into the matter era and after the transition 
the radiation disappeared. Under this assumption one gets the simple expressions t^^"^ 

^The evolution equation for a{t) may be transformed into an equation for a{z) and then for f[R(z)] 
provided z{t) is a monotonic function. As discussed above this means that one takes into account only 
this class of solutions. It is unclear as to what extent this restriction affects the final outcome. 
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and t^/^ for a{t) in these eras respectively; this is the basic assumption in 



It should be stressed that in general it is impossible to reconstruct a differential 
equation from its approximate solution. A given function is an approximate solution to 
a whole class of equations and the correct equation is in no way distinguished within this 
class: the Einstein equations correspond to the saddle points rather than to attractors 
of the dynamical system of ref. Metric NLG theories are so conceptually involved 
that this (rather obvious) fact goes overlooked. As an example we take a simpler theory. 
Recall that it is possible to reconstruct the Newtonian theory (Laplace equation) from 
the exact solution, the Newton's potential. Now one tries to reconstruct it from an 
approximate solution. Suppose an astronomer observes a planet of mass m moving on 
a very elongated elliptic orbit with the eccentricity e = l — £, £:<^1 and e positive. He 
assumes that the gravitational force comes from a central potential U{r) and tries to 
determine it from the planet's motion. To this end he employs the energy integral 



where J is the conserved angular momentum. The astronomer determines from obser- 
vations the positions r at different times and fits the data by an analytic function r(t), 
inverts it and expresses the derivative f as a function F{r). He makes his observations 
in two regimes: when the planet is close to perihelion and to aphelion. From the exact 
solution of the Kepler problem we may compute that near the perihelion the dependence 
t[r) for this elongated orbit is approximately (by keeping only the dominant term in e) 



for U = and this function should be used by the astronomer as an analytic fit to 
his data set. When the planet moves near the aphelion the approximate solution t{r) 
which also fits the other observational data is 



Both the functions are easily invertible and each is used to compute r = F{r). The 
quantities a (the semimajor axis), e, J and the parameter p are determined from the 
observations (assuming ellipticity of the orbit), then en = ^ and upon inserting F(r) 
into the energy integral the astronomer finally finds that for r close to ry^\„ 



2 zmr'^ 





U{r) 



J2 



J2 



+ E 



mpr 



®In conformity with classical mechnics we denote it by U instead of 0. 
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while close to the aphelion 



U{r) 



J2 



(2a - r) 



J2 



+ E. 



2mr2 



None of these expressions is an approximation to the Newton's potential. An ultimate 
conclusion from this example and the two works mentioned above is that reconstructing 
the underlying Lagrangian from known Robertson-Walker solutions is impossible for 
fundamental reasons. Relativistic metric theories of gravity are too intricate. It is in 
order here to remind the Einstein's view that a new physical theory is never formulated 
by induction from a pile of empirical data. 

All that above does not imply that alternative theories of gravity should not be 
applied in cosmology and in particular should not be used to account for the dark en- 
ergy. Today the situation in gravitational physics is exceptional as compared to other 
branches of physics: the well established and confirmed theory, general relativity, seems 
to be just a point in the "space" of all existing and conceivable theories of gravitational 
interactions and its nearest neighbourhood is densely populated by its alternatives, the 
metric NLG theories. According to Cantor's theorem the cardinal number of the set 
of metric NLG theories is greater than the cardinality of the continuum. The very 
existence of these theories entitles one to apply them to describe effects which are grav- 
itational or may be interpreted as such. On the other hand the wealth of these theories 
makes necessary, before making any applications of one chosen from this huge set, to 
investigate two problems: i) to determine all possible interrelations between them and 
their relationships to general relativity, ii) whether a given gravity theory satisfies all 
the well grounded general rules of classical field theory and has acceptable properties. 
These include: 

— determination of particle contents (spectrum), 

— existence of a stable maximally symmetric ground state, 

— form of interactions with ordinary matter, i.e. which quantities are measurable. 
A natural and obvious postulate is to require that a given gravity theory should have 
a Newtonian limit. It however creates conceptual problems. It is not sufficient to for- 
mally get in the linear approximation the Poisson equation for a scalar potential [30] 
Af/ = AirGp + A or another equation appearing in the Newtonian gravity (such as the 
stellar hydrostatic equilibrium one [SI])- Newtonian gravity is well defined as a small 
static perturbation of Minkowski spacetime and the postulate may be unambiguously 
formulated for those Lagrangians which admit the fiat spacetime as a solution. In this 
class of theories the postulate may work as an effective discriminating condition. Yet 
the NLG theories are most interesting and attractive for cosmologists if their ground 
state solution is curved being de Sitter or anti-de Sitter space while the fiat spacetime 
is excluded. In this case the very concept of Newtonian gravity in these spacetimes is 
controversial and there are confiicting results on the weak-field limit for different La- 
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grangians. As far as we know the problem of the Newtonian hmit is not convincingly 
solved even in general relativity in the presence of the cosmological constant (a pre- 
cise definition of the limit is lacking). Following cosmologists we deal mainly with the 
theories which do not admit Minkowski space as a solution and at the present level of 
knowledge the Newtonian limit criterion cannot be applied to them. 

Investigations of the second problem, i.e., to what extent the metric NLG theories 
are viable from the viewpoint of classical field theory are the main subject of this paper 
and the subsequent Paper II. In section 2 we count the number of degrees of freedom 
and briefly discuss the particle spectrum of a general NLG theory. The latter subject 
requires replacing the fourth-order Lagrange equations by dynamically equivalent second 
order ones for the resulting gravitational triplet of fields; this aim is achieved by the 
powerful method of Legendre transformations. This decomposition is crucial for all 
investigations of the theory. The triplet may be described in infinitely many different 
frames and the two most important ones and the problem of coupling ordinary matter 
to the gravity fields are presented in section 3. The controversial problem of which 
frame is physical does not affect the criteria of viability of various gravity theories. The 
most fundamental criterion is the existence of a stable maximally symmetric ground 
state solution and it is investigated in Paper II. Conclusions and further considerations 
concerning the possibility of recovering the underlying L — f{R) Lagrangian from a 
given (cosmological) solution are contained in section 4. 

2 Particle spectrum 

We consider in this section the general metric NLG theory based on an arbitrary La- 
grangian L = f{gfj,u,Rai3^iu)- The theory is metric in the sense that a nondegenerate 
tensor field g^i, with Lorentzian signature is the only independent dynamical quantity. 
One may also investigate a metric-affine theory ("the Palatini method") with the same 
Lagrangian wherein one takes independent variations of L with respect to the metric and 
a symmetric connection. As is well known, for any L different from R the two theories 
diverge. Some authors claim that the metric-affine approach is more natural since the 
Lagrange equations are of second order while in the purely metric theory they are of 
fourth order. However the metric NLG theory is not inherently a higher derivative one. 
The tensor g^^i, appearing in the Lagrangian actually is a kind of unifying field mixing 
various particles (fields) with different spins and masses. To find out a physical inter- 
pretation of the field it is necessary to decompose it in a multiplct of these fields. Then 
equations of motion for the separate fields are of second order and display a physical 
content of the theory better than the original fourth-order ones. This is why referring to 
NLG theories as "higher derivative ones" is misleading. In this respect the metric-affine 
approach is not advantageous over the purely metric formalism and we prefer the latter 
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as conceptually simpler. 



To avoid any confusion and for sake of completeness we begin with determining the 
degrees of freedom of the general metric NLG theory though it may be found in the 
literature. Counting the degrees of freedom (d.o.f.) in this case is far from trivial. The 
quadratic theory, L = R + + R^.^R^''' , is known to have eight d.o.f. [D E21 El El] • For 
Lagrangians with arbitrary dependence on the Ricci tensor and applying a perturbative 
approach to Lagrangians depending on the Weyl tensor the d.o.f. were first counted 
by Hindawi et al. [35] by using a second order version of the theory. For most general 
Lagrangians it is possible to determine the maximal number of d.o.f. in the initial 
fourth-order formulation. It is well known that the pure gravitational field in general 
relativity, R^^ = 0, has two d.o.f. [30] and we shall count them in the same way for 
L = f{gfj.u,Rai3fiu)- The number of d.o.f. for a given system is defined as a half of the 
number of arbitrary functions needed to uniquely specify the initial data for the Cauchy 
problem for the equations of motion of the system. The equations of motion following 
from this Lagrangian form a system of ten tensor PDE of fourth order E"^ = for the 
unifying field g^^. Let a spacelike hypersurface S be chosen as an initial data surface 
for the equations. The theory is generally covariant (diffeomorphism invariant) and one 
can freely choose a coordinate system in the spacetime and the most convenient one is 
the comoving system (normal Gauss coordinates), qqq = —1, goi = 0, such that S has an 
equation t = 0. There are six unknown functions gik and the initial Cauchy data consist 
of values of gik, dgik/dt, S^gik/dt^ and d^gik/dt^ on S; these are 24 functions of three 
coordinates x\ The data are subject to a number of constraints. First, the coordinates 
on S may be freely changed and this gives a freedom of choice of 3 functions of these 
variables. Secondly, the trace of the extrinsic curvature of S may be given any value, 
i.e. one function is arbitrary. Finally, the general covariance of the theory implies, in 
the same way as in general relativity, that some of the field equations are constraints. 
In fact, the invariance of the action 



under an infinitesimal coordinate transformation gives rise to a strong Noether conser- 
vation law ("a generalized Bianchi identity") = (see e.g. appendix A in arXiv 
version of [37]). Explicitly it reads 



The last three terms contain at most fourth time derivatives of g^^, and the identity 
implies that OqE^ cannot involve fifth time derivatives. Thus involve at most third 
time derivatives and E'^ = are not propagation equations but form four constraints on 
the initial data. Together the number of independent Cauchy data is diminished to 16 
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arbitrary functions and thus a general metric NLG theory has eight d.o.f. 



For a restricted NLG theory, L = f{R), the number of d.o.f. is less than 8 what 
means that the equations = generate additional constraints on the Cauchy data. It 
is not easy to determine all the constraints from these equations and one should instead 
apply a second order formulation of the theory. To this end it is adequate to view the 
L = f{R) theories in a wider context of as large class of of NLG theories as possible. 
We therefore consider for the time being Lagrangians L = f\g^u, Rap), i-e. with no de- 
pendence on the Weyl tensor. 

An adequate mathematical tool for this purpose is provided by a specific Legendre map 
[38l [39| HO] . For these NLG theories the method is as general and powerful as Legendre 
maps transforming the Lagrangian formalism into the Hamiltonian formalism in clas- 
sical mechanics and classical field theory. Yet the method is not currently used in a 
systematic way and most papers on applications in cosmology have employed various 
ad hoc tricks to transform from the Jordan frame to the Einstein frame. The tricks in 
most cases give results equivalent to the Legendre transformation, however do not allow 
to fully display the structure and features of the theory. Here we give a brief summary 
of investigations of particle spectrum contained in [H] while the general formalism is 
described in [88| [89t [¥)] . 

The Jordan frame consists of only one dynamical variable, the tensor field g^^, JF = 
{9pLi>}i which plays both the role of a metric tensor on a spacetime M and a kind of 
unifying gravitational field being a composition of some fields having definite spins and 
masses. Pure gravity is then described by a multiplet of the fields having together at 
most eight d.o.f., the metric is a geometric component of the multiplet. The unifying 
field may be decomposed into the component physical fields in two ways. The first 
method assumes that g^^ is the spacetime metric and one separates from it, by means 
of a Legendre map, the additional degrees of freedom, i.e. the other components of the 
multiplet. The Ricci tensor is decomposed into its irreducible parts, the trace R and 
the traceless tensor Sfj,y = R^u — \g^uR, then one defines a scalar and a tensor canonical 
momentum conjugate to the "velocity" R^v by Legendre transformations 

Together with the metric the two fields form a gravitational triplet, named the Helmholtz- 
Jordan frame, HJF = {gp.viX-,'^^^}- Equations of motion for the triplet follow from a 
Helmholtz Lagrangian Lh which is dynamically equivalent to the initial L = f{gij,u, Rap)- 
Here one meets a technical obstacle: to get an explicit form of Lh one must express R 
and S^y in terms of x and n^^" , i.e. to solve the defining equations to find R = r(x, vr^'^) 
and S^u = 7r°^). For a general L = f{g^u,Rap) this requires solving nonlinear 

matrix equations. (No doubt, the power of Hamiltonian formalism in physics stems 



13 



from the fact that physically relevant Lagrangians are quadratic in "velocities".) This 
is why only Lagrangians that are quadratic in the Ricci tensor have been investigated 
in detail. However the formalism in principle works for any L. Furthermore, as we shall 
see, in the Einstein frame the specific dependence of L on R^^ only affects interaction 
terms while the general structure of the theory remains unaffected, therefore although 
a general formalism (in HJF) has been developed for any function /, in practice one 
applies only those / for which the Legendre transformations may be effectively inverted. 
For any / the Helmholtz Lagrangian reads 



where if is a Hamiltonian. This Lagrangian is linear in R^^u, what implies that Lagrange 
equations for ^f^,^ take form of Einstein field equations, G^j/ = 8nGTfj,ty{Rai3, x, vr°^). The 
RHS of these equations is by definition an energy-momentum tensor for the two fields 
and it depends linearly on the Ricci tensor and on the first and second derivatives of 
X and TT^*^. There are no kinetic terms for x aiid tt^'^ in Lh and propagation equations 
for the fields are derived in a rather intricate way from T^^, these are hyperbolic second 
order ones. Both the fields are subject to one algebraic and four first order differential 
constraints and in consequence they carry together six degrees of freedom. In summary, 
the particle spectrum of the theory exhibited in HJF consists of: a massless spin-2 field 
(gravitoE0, spin two and 2 d.o.f.), a massive spin-2 field (5 d.o.f.) and a massive scalar 
field. This outcome (with the same values of the masses for the two fields) was first found 
in the linear approximation [U HS] and then by various methods in the exact theory for 
a quadratic Lagrangian R + R"^ + R^^R^'^ [SI HSj. HJF is not uniquely determined: 
since the Helmholtz Lagrangian in (1) is not in a canonical form, various redefinitions 
of X and TT^'^ are admissible. 

Introducing the field n'^'^ makes sense if its definition may be (at least in principle) 
inverted to yield S^u = s^i^ixi tt"^), otherwise vr^'' = 0. The field exists if the Lagrangian 
depends on R^^ in a nontrivial way, i.e. the Hessian 



This condition does not hold for L = f{R)- Then vr^'^ vanishes and gravity is described 
by a doublet HJF = {gfj,u,x} carrying 2+1 d.o.f. and the scalar cannot be massless. 
Thus for this class of Lagrangians determining the particle spectrum is very simple and 

^We use this traditional name for the metric field satisfying R^^, = 0. It is worth noting, however, 
that the relationship between the hypothetical quantum of the gravitational field and the classical field 
(described by general relativity) is different from that between the photon and the classical Maxwell 
field, see 



Lh = R + xR + ^^"S^, - H{x, Tin 



(1) 




14 



straightforward. 



The other approach to constructing a second order formahsm is more sophisticated. 
Here one assumes that g^i, is merely a unifying field for gravitation and plays a role of a 
spacetime metric in a purely formal way — in the sense that the Ricci tensor appearing in 
the Lagrangian is made up of it. One introduces a new metric as a canonical momentum 
conjugate to the full Ricci tensor via a Legendre map as 
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df 
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a/3, 
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df 



dR 
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here g = det{g^u). g^^ may be viewed as a metric tensor providing that detidf/dRap) 7^ 
0. Clearly for f{R) = R one gets g'^" = g^^ and for arbitrary f{R) the new metric is 
conformally related to the old one, g^y = f'{R) g^u, where is the matrix inverse to 
g^'^. In this case the Legendre transformation is degenerate since it cannot be inverted. 
In general the transformation is a map of a (Lorentzian) metric manifold {M^g^i,) into 
another one, {M,g^j_y). If the transformation (2) is invertible, i.e. the Hessian for / does 



not vanish, it may be solved to give R 
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As in Helmholtz- Jordan frame 



one constructs a Helmholtz Lagrangian which now takes a generic form 
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f{gf,u,ra(3ig,g)), 



(3) 



where V is the covariant derivative with respect to It is worth stressing that 

in Einstein frame, EF = {gfiu,gai3}, precisely the Einstein-Hilbert Lagrangian for the 
spacetime metric is recovered, giving rise to Einstein field equations G^uijj) = 
8TTGTf^i^{g, g) with the tensor field g^^, acting as a "matter" source for the metric. The 
tensor T^j, is the variational energy-momentum tensor for ^f^jy defined in the standard 
way and contains second derivatives V/S^^gap but no curvatur^. Unlike Lh in HJF, the 
Lagrangian has a canonical form for gaf3, i.e. is a sum of a kinetic term and a potential 
part. The kinetic term, KCVg), is a quadratic polynomial in first derivatives V^gap, 
as is usual in classical field theory and what is really remarkable, it is universal, i.e. is 
independent of the form of the function / |3H1 EH]- The only reminiscence of the original 
L in JF is contained in the potential part of Lh'- explicitly via fig^iv, Tap) and implicitly 
via Tap, i.e. in interaction term^. Hence in EF one recovers just general relativity with 



^In general the energy-momentum tensor involves second derivatives of matter variables, the gauge 
fields, the minimally coupled scalar field and perfect fluids belong to few exceptions. 

^Clearly this is not little. In quantum mechanics every state vector satisfies the Schrodinger equation 
and the whole variety of quantum systems is encompassed in interaction terms in the Hamiltonian. 
Here something analogous occurs. We stress this point since in the fourth-order formulation of an 
NLG theory in JF an impression arises that the theory is more different from general relativity than it 
indeed is. 
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a source which may be interpreted either as a nongeometric component of gravity or 
merely as a (quite exotic) matter field described by classical field theory. In this sense 
Einstein general relativity is a universal Hamiltonian image (under the Legendre map) 
of any L = f{g^i,, Rap) gravity theory. In other terms, general relativity is an isolated 
point in the space of all gravity theories: its closest neighbourhood, consisting of the 
metric NLG theories, can be mapped onto it and thus is not different from GR. It is 
also clear that in practice there is no need in studying Lagrangians more complicated 
than quadratic in R and R^u- 

The second order Lagrange equations for g^^i, in EF are subject to four differential 
constraints (following from Bianchi identities VuG'^^{g) = 0) which allow one to elim- 
inate four of ten components of the field. This shows that it carries six d.o.f. and is 
actually a mixture of two different physical fields. The next step is thus to decompose 
it into components with definite spins. Then one again gets a scalar and a spin-2 field. 
For practical purposes it is convenient to eliminate the scalar from the outset by an 
appropriate choice of the Lagrangian in JF [H]. The field ip^i, arising in this way from 
g^i, carries spin two (and five d.o.f.), has the same mass as that computed in HJF for n^'^ 
and is nonlinear (it is well known that any linear spin-2 field is inconsistent in general 
relativity [IS]). It is straightforward to show in EF that ip^u is necessarily a ghost field 
(a "poltergeist") [U HSl HH US, HT] , while it is rather difficult to estabhsh this feature for 
TT^'^ in HJF. And to avoid any misunderstanding we stress that the ghost-like behaviour 
of the spin-2 field is inevitable: it appears in any consistent theory of gravitationally 
interacting spin-2 fields |1H] and in particular is a feature of any L = f{g^u, Rap) gravity. 

All these gravity theories are similar and they differ only in the interaction terms 
in (3) and in masses of the spin-2 and scalar components of the gravitational triplet. 
All the L = f{R) theories are reduced in Einstein frame to general relativity plus a 
massive minimally coupled field with a self-interaction potential determined by / [H7] . 
At first sight it seems that Lagrangians generating tachyonic masses of the fields should 
be excluded as untenable. This is the case when the ground state solution is Minkowski 
space and other fields behave as small excitations in this spacetime. However, if the 
ground state is anti-de Sitter, also the scalar field with a tachyon mass is allowed if its 
modulus is not too large in comparison to the cosmological constant, see Paper II. 

The case of most general Lagrangians, L = f{gfj,u, Rai3,Cfiuai3), explicitly depending 
on the Weyl tensor, is more subtle. The unifying field g^^ carries 8 d.o.f. and one 
conjectures that it can be decomposed in the same triplet of the graviton, a spin-2 
field and a scalar, however a proof for an arbitrary L is missing. The conjecture was 
proved only in a perturbative analysis (and in four dimensions): one expands a generic 
Lagrangian about a ground state solution of the theory (Minkowski, de Sitter or anti-de 
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Sitter) up to terms quadratic in Riemann tensor and arrives at HH] 

L = const + R + aR'^ + bR^^R'"' + cGB, 

where GB, the Gauss-Bonnet term, is a topological invariant and in d = 4 may be 
discarded as a divergence. Thus the most general NLG theory perturbatively reduces 
to that without the Weyl tensor and has the same particle spectrum. 

We stress again that in both the frames any NLG theory is reduced to general 
relativity plus some exotic source. In EF the general-relativistic form of the theory 
is obvious, in HJF it is less conspicuous due to the specific form of the Helmholtz 
Lagrangian (1). From the physical viewpoint the spin-two and spin-0 fields may be 
viewed either as the components of the gravitational triplet or just as a kind of matter. 
It is not quite clear whether the difference between the two interpretations is empirical: 
whether there is a "gedankenexperiment" allowing to differentiate one from the other. 
As in Brans-Dicke theory it is necessary to assume that the two fields do not couple 
to any other matter in the sense that in a relevant Lagrangian there are no interaction 
terms of the two fields with the particles of the standard model. And as in Brans-Dicke 
theory "one way" interactions are admissible: ordinary matter may act as a source in 
equations of motion for x and vr^'^. The two nongeo metric components of gravity only 
interact gravitationally in the sense that there is interaction between them and, first of 
all, they act as a source of the spacetime metric in Einstein field equations. It is often 
assumed (and clearly there is no proof) that the dark energy signals its existence solely 
by its influence on the cosmic evolution. Whether it should be regarded as a form of 
matter or as a component of gravity is presently a matter of convention. Mathematically 
the issue is irrelevant and Einstein field equations arising from (1) and (3) should be 
studied by applying all the methods developed to this aim in general relativity. In 
particular, the energy-momentum tensors for the spin-0 and spin-2 fields appearing in 
these equations ought to satisfy the conditions usually imposed on matter in general 
relativity. Clearly it may be claimed that the two components of gravity are specific in 
the sense that they need not satisfy the standard conditions, e.g. the scalar field may not 
be subject to energy conditions. It might be so, however the price would then be high: 
the most significant results (or their appropriate analogues) found in general relativity 
would be inaccessible in a given NLG theory. In this work we assume that the spin-0 
gravity is not exceptional in this respect. 

3 Frames and initial conditions 

Once an NLG theory is expressed in HJF one may perform arbitrary Legendre trans- 
formations ("canonical ones") and changes of variables (field redefinitions), thus the 
theory may be formulated in infinitely many various frames. Clearly both HJF and 
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EF are priviledged by their construction exhibiting the physical content of the theory. 
All the frames are mathematically equivalent provided the transformations are at least 
locally invertible (proving a global invertibility is a hard task) and nonsingular. Equiv- 
alence means that the space of solutions in one frame is in a one-to-one correspondence 
to the space of solutions in another frame. The corresponding solutions are different 
and physical quantities made up of them are different, most notably energy is very sen- 
sitive to various transformations. Thus dynamical equivalence of frames implies their 
physical inequivalence. As long as the theory is closed, i.e. all that exists is contained 
in a Lagrangian of the theory (in the case discussed here "everything that exists" is 
pure gravity and there is no matter), this inequivalence is irrelevant as undetectable. 
All the frames are equally physical. For example, spacetime intervals between a given 
pair of events are different in distinct frames and the differences cannot be measured 
without external rods and clocks. And to measure energy of the gravitational triplet 
one needs an external device which is not included in the Lagrangian, yet in a closed 
theory neither an external observer nor external device does exist. In the same way 
spacetime intervals between a given pair of events are different in different frames and 
their differences cannot be measured without external rods and clocks. In this sense 
all the theories of physics are open: the observer and his equipment is not described 
by a tested theor}o. To make an NLG theory open it is necessary to couple it to or- 
dinary matter and predict then some effects which may be observed by an external agent. 

Coupling of matter to gravity should proceed in the same way as in general relativity 
where, however, no ambiguity appears since there is only one frame. For concreteness 
we consider now the restricted NLG theories, L = f{R), since we will be dealing with 
them in Paper II. Clearly the coupling of matter to gravity is the same in any NLG 
theory. One takes pure gravity and chooses a frame consisting of a tensor 7^,^ regarded 
as a spacetime metric and a scalar 0; these quantities are some functions of the variables 
g^y and x forming HJF. The corresponding Lagrangian may have almost arbitrary form 
Lg[R['y),(j)) [m [37]. For a given kind of matter its Lagrangian is constructed in 
special relativity and gets some form Lm{ri^^,'^,d'^). By definition, 7^1. is the metric 
of the physical spacetime and any matter minimally couples to it, hence the matter 
Lagrangian for becomes Lm(7^,^, \E', V\I'), where V is the Levi-Civita connection for 
7^j,. The scalar gravity (p does not couple to the matter and the total Lagrangian is 
just Lt = Lg + Lm- The chosen frame is the physical frame due to the minimal cou- 
pling. In all other frames where the transformed tensor playing the role of a metric is 
different from 7^j,, the matter is nonminimally coupled to it and a coupling to may 
appear. Dynamical equivalence of various frames remains preserved in presence of any 
matter while these frames should be regarded as unphysical since experimental devices 

^''This is not so trivial as it may seem. There are some tendencies in quantum gravity to regard it 
as a closed theory. 
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measure quantities made up of variables of the physical frame. For example, optical 
observations disclose that the light of distant galaxies is redshifted, what is interpreted 
as that expanding Robertson-Walker spacetime forms the physical frame, while one 
may make all computations in a conformally related frame where the spacetime is flat. 
In this case outcomes of the computations must be transformed back to the physical 
frame, the R-W spacetime, if they are to be confronted with observations. There is 
nothing new in this, the same is always done in classical mechanics: first one determines 
(experimentally) physical positions and momenta of a given system, makes a canonical 
transformation mixing these quantities to new variables in which the Hamilton equa- 
tions are easiest solvable and finally makes the inverse transformation to express a given 
solution in physical variables. In classical mechanics this is obvious, in a gravity theory 
it is not. 



Assume, as most authors applying NLG theories to cosmology actually do, that JF 
is physical, then the total Lagrangian is Lt = f{R) + Lm{g^u, ^, V\l/) with V being now 
the metric connection for g^^. The Helmholtz Lagrangian in HJF reads 

Lh = p[R{g) - r{p)] + /(r(p)) + L^{g,,, M/, VvP), (4) 

from now on the spin-0 gravity in HJF is denoted p and defined as p = ^ in conformity 
with [HHl ST]. The Lagrange equations are then R{g) = r{p), 

Gf,u{g) = 0^,{g,p) + ^t^uig, ^) (5) 

and ^ = 0. (6) 

Here r(p) is a (possibly unique) solution of the equation = p and 9f^u is an effective 
energy-momentum tensor for p [37] while t^y{g, \E') is the standard variational energy- 
momentum tensor for matter derived from with the aid of g^j.^. A propagation 
equation for the scalar is derived, as previously for x, from the tensor 9^^ by taking the 
trace of (5) and it reads 

Up - ^/(r(p)) + ^pr{p) = -^gf^-^t^^ig, ^). (7) 

It is easy to see that the scalar does not have its own Lagrangian and eq. (7) must 
be derived in this roundabout way. The metric field has two sources, the matter and 
spin-0 gravity. The equation of motion for any matter, (6), is independent of p, yet the 
matter forms a source term g^^t^j^u for the scalar gravity. In this sense the field p does 
not directly affect motions of matter and its effects are confined to affecting the metric 
via eq. (5). 
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The initial and/or boundary conditions for both g^^i, and p are related to or even 
determined by a matter distribution. For example, in the solar system the matter dis- 
tribution is dominated by the sun and both the fields are static spherically symmetric. 
Inside the sun the strength of the material source for g^u and p is comparable and both 
the fields are there relatively strong and regular at the centre. Outside the sun the fields 
are fading towards spatial infinity. Notice that the scalar generates spherically symmet- 
ric corrections to Schwarzschild solution. Search for these corrections has resulted in the 
fact that the Lagrangian R + l/R is ruled out by measurements of the PPN parameter 
7 in the solar system |T2] and provides very stringent bounds on the size of possible 
corrections to the Einstein-Hilbert Lagrangian L = R [191 E [IE]. This is an indication 
that Jordan frame is unlikely to be the physical frame. 

If instead, Einstein frame is regarded as physical, the Helmholtz Lagrangian takes a 
form well known from general relativity, 

Lh = R{g) - - 2V(0) + L^{g, v[/, Vvl/), (8) 

where = -y/llnp, is a potential determined by / and matter is minimally coupled 
to g^^. The Lagrangians in (4) and (8) have the same dependence on the spacetime 
metric, g^^p and g^p, respectively. As a consequence the matter energy-momentum tensor 
in EF is just tf^^{g, \E'). The field equations are now directly derived as the variational 
ones, 

G^M=T^A9,<P)+t^A9,'^) (9) 

a,. IK ,0) 

and = 0. (11) 

Obviously one obtains eq. (11) from (6) by replacing g^i, with Now the scalar 
gravity is completely decoupled from matter and solely interacts with the metric field. 
This implies that initial and boundary conditions for are independent of matter dis- 
tribution. The metric has two independent sources and its symmetries, boundary and 
initial conditions are determined by both or by the source that dominates. For example, 
in the solar system the field need not be spherical. From the spherical symmetry of 
the local spacetime one infers that the matter (the sun) is dominant on the RHS of eq. 
(9) and any small deviations from Schwarzschild metric should be nonspherical. The 
scalar field has no local (matter or any other) sources and fills the entire universe and 
approximately is described by one global solution. The solution is approximate since 
the scalar interacts and is affected by the spacetime metric and the latter is affected 
by local matter inhomogeneities. This global solution which is realized in our universe 
was determined by some initial and boundary conditions near the Big Bang. From the 
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observations showing that the cosmic space is homogeneous and isotropic on large scales 
it is inferred that either the (f) field is homogeneous (and time dependent) throughout 
the spacetime since the Big Bang or it is inhomogeneous and undetectably (even using 
cosmological data) weak. Clearly the first possibility is more attractive as it may a 
priori account for the dark energy. In any case a solution for the scalar has been chosen 
once for the universe and unlike the metric field it cannot be locally varied according 
to local (arbitrary) conditions. In other words, even in vicinity of a black hole, where 
the spacetime is extremely distinct from the Robertson-Walker one, the field only to 
some extent deviates from its overall cosmological solution. We emphasize that if some 
nonspherical deviations from the Schwarzschild metric are detected in the solar system 
(after subtracting all effects of planets, the Kuiper belt, the Oort cloud etc.), they may 
be accounted for by the scalar component of gravity in Einstein frame. 

Both the frames are experimentally distinct. Which of them (if any) is physical? In 
other terms, which metric is minimally coupled to matter? The problem arises in any 
theory in which various frames appear. For example, in string theory in the low energy 
field-theory limit of string action one may use either the string (Jordan) frame in which 
the stringy matter is minimally coupled to the metric while the dilaton field is nonmin- 
imally coupled to it or transform to the conformally related Einstein frame where the 
dilaton is minimally coupled to the new metric and has the canonical kinetic term. The 
two frames are usually considered as completely equivalent for describing the physics of 
the massless modes of the string. This is particularly noticeable in the pre-Big Bang 
infiationary string cosmology. A superinfiationary solution in string frame becomes an 
accelerated contraction in Einstein frame and vice versa. This drastic difference in be- 
haviour of the cosmic scale factor in both the frames is irrelevant for string cosmology. 
In fact, the number of strings per unit of string volume is decreasing in time during the 
pre-Big Bang inflation in both frames and the temperature of the string gas grows in 
comparison to the temperature of the photon gas in both frames [HD]- Whether a{t) is 
expanding or contracting, the horizon/flatness problem of cosmology may be solved in 
each frame. 

This frame independence for the physical effects of inflationary solutions in string cos- 
mology, showing invariance of physics under local fleld redeflnitions, has been found, 
however, only for a limited number of observables among those which can be constructed 
in this theory. Other observables are frame dependent. Even in the early universe space- 
time intervals are measurable, at least in principle if not in practice, and these quantities 
depend on whether a{t) grows or decreases. The very fact that we observe that both 
the optical spectra of distant galaxies and the cosmic microwave background radiation 
are redshifted, indicates that cosmology is not frame independent. The pre-Big Bang 
era of string cosmology is so far from us and exotic that very few effects may be ob- 
served today and these turn out to be frame invariant. Yet it is doubtful if string theory 
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as such is frame (i.e. conformally) invariant!"]. This would mean that the theory is closed. 

In metric NLG theories which are closer to experimental ph ysi cs and observational 
astronomy, the fact that most observables are frame dependenlo, i.e. only one frame 
is physical (is unique up to trivial field redefinitions), is hardly arguable (see however 
references in [37]). The ultimate decision of which frame is physical will be given by 
experiment, but it should not be expected very soon. Before making any application 
of these theories it is necessary to establish by some theoretical arguments or by mere 
assumption which frame is measurable. More than ten years ago we studied in detail 
the full network of relationships between restricted NLG theories, scalar-tensor gravity 
theories and general relativity and effects of introducing matter in various frames and 
on this basis we gave arguments in favour of Einstein frame After our work there 
was much discussion in the literature on the subject [H2]. Since the discovery that a 
modified gravity might replace the dark energy most authors have preferred the Jordan 
frame as physical. As it is motivated in a paper, "if one wants to consider modifications 
of gravity like scalar-tensor theory or metric f{R) gravity, the Jordan frame should be 
assumed to be the physical one" since if Einstein frame is physical "the resulting theory 
will be no different from general relativity" . These authors seem to be unaware that also 
in Helmholtz- Jordan frame the field equations for the metric tensor are Einstein ones. 
The assumption that by choosing Jordan frame as physical (and minimally coupling all 
matter in this frame) one gets a theory which is essentially more different from general 
relativity than in the other case, is a mere illusion. 

We do not wish to enter the debate again. In our opinion all relevant arguments in 
favour of JF and EF have already been expounded. We only point out that in most 
cases the Einstein frame is computationally advantageous (and this is explicitly or im- 
plicitly acknowledged by most authors who implicitly or explicitly assume Jordan frame 
as physical). For example, the Zeroth Law and the Second Law of black hole ther- 
modynamics for a polynomial f{R) have been proved only in Einstein frame [53]. On 
the other hand we keep open mind for the possibility that the real world might be not 
so simple as we expect. It might be so that the real physics is in Jordan frame while 
Einstein frame remains advantageous both in solving equations of motion and, what is 
most important, in proving general features of a given theory. The issue will be resolved 
by experiment. 

The reader may find these comments on physical interpretation of various frames as 

^"'^In string theory, which is still far from making concrete physical predictions, it is difficult to establish 
what is measurable. Some authors claiming complete frame invariance of the theory seem to confuse 
the dynamical equivalence of various frames, which is indisputable, with physical equivalence. 

"'^^Nevertheless in Jordan and Einstein frames many physical quantities are the same, e.g. for black 
holes all the thermodynamical variables do not alter under a suitable Legendre map [51] . 
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unnecessary. We explain that we make these comments exphcitly and in some length 
since for many researchers in cosmology these things are far from being clear. The level 
of confusion may be best seen from the introduction to the work [54j where the authors 
express their surprise that the following procedure: i) transform the equations of motion 
from Jordan frame (which is assumed to be physical) to Einstein frame, ii) solve the 
transformed equations and choose these solutions which seem to be physically signif- 
icant and iii) transform back the chosen solutions to JF, does not provide physically 
meaningful results. They even quote a paper where it is claimed that "passing from 
one frame to the other can change the stability of the solution" — a statement which 
for a mathematically oriented reader clearly means that the transformation is singular. 
The authors of |5l] summarize their doubts in the question "is it correct to obtain a 
solution in a frame and then interpret it in another frame?" . To leave no space for any 
ambiguity we respond: yes it is, in the same sense as it is done in classical mechanics 
where we deal with canonical transformations. 

We emphasize that for the purposes of the present work the issue of which frame is 
physical, is irrelevant. We claim that the order of investigations should be as follows. 
First one chooses (on some basis) one or a class of NLG theories. Second, one veri- 
fies if the chosen theory is viable from the viewpoint of classical field theory. Third, 
a physical frame is assumed. Then equations of motion are solved either in JF or in 
HJF, EF or some other frame, depending on computational facilities (usually EF is most 
convenient). Finally the solutions should be transformed back to the physical frame (if 
were found in another one) to construct physical observables which will be confronted 
with observations. We stress that the second step cannot be passed over otherwise there 
is danger that the work will be wasted. In fact, in some papers long and nontrivial 
computations have been performed in a framework of a theory which is definitely un- 
tenable. Viability criteria are independent of the choice of the physical frame. This 
frame independence does not mean that whether a given theory is physically viable or 
not can be established in any frame. On the contrary, all the methods developed up 
to now to investigate the viability do work in Einstein frame and usually (besides one 
case mentioned in Paper II) do not work in most other frames. In this sense Einstein 
frame is mathematically distinguished. Yet these methods neither prove that this frame 
is physical nor assume it. We shall use Einstein frame for checking viability various 
L = f{R) theories. Therefore the first, basic assumption or criterion a theory should 
pass is the existence of Einstein frame: the Legendre transformation from JF to EF 
must be regular (in a neighbourhood of a candidate ground state solution). Once a the- 
ory meets the criteria and is regarded tenable one may choose the physical frame. We 
shall see that there are infinitely many viable theories and untenable ones are equally 
numerous. Obviously of all viable restricted NLG theories only one (if any) is correct, 
i.e. will be fully confirmed by experiment and observations. By the time this occurs one 
may view any viable theory as a candidate for describing gravitational interactions in 
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the nature. 



Out of all possible viability criteria for a classical field theory the most appropriate 
(and effectively applicable) one in the case of metric NLG theories is the existence of a 
stable maximally symmetric ground state. Paper II is devoted to studying and applying 
this criterion. 

4 Conclusions 

It is reasonable to conjecture that the accelerated expansion of the universe is not driven 
by an extremely exotic and unknown to physics kind of matter with negative pressure 
but is rather due to some modification of gravitational interactions. However in search 
for a modified gravity theory great caution is necessary. Gravitational physics is ex- 
ceptional among all branches of physics in that there is a great variety of competing 
theories, all of which are some variations of Einstein's general relativity. Modifications 
may go in all possible directions while in most cases gravity theories applied to cosmol- 
ogy differ from general relativity only in one axiom: the form of the field equations. The 
assumption that a gravitational Lagrangian is an arbitrary function L = f{gfJ,u,Ra|3^lu) 
gives rise to infinity of theories and choosing the correct one is a hard task. It is almost 
invariably attempted to make this choice employing the cosmic scale factor a{t) in the 
spatially fiat Robertson- Walker spacetime. This approach is actually hopeless even in 
the framework of restricted metric gravity theories with L = f{R). It has been shown 
[28] that if one knows the exact analytic form of a{t) and a number of initial condi- 
tions, then the function f{R) is a unique solution of a linear third order ODE. This 
method does not work in practice since the equation is intractably complicated. The 
fundamental cause that the method cannot work is that the real physical spacetime is 
not R-W one. The cosmic matter distribution becomes homogeneous and isotropic only 
asymptotically at large scales. If one were ingenious enough to solve the differential 
equation for f{R) employing that form of a{t) which best fits all the astronomical data, 
the resulting Lagrangian would be rather different from the true one. The standard 
Friedmannian cosmology based on general relativity is the best fit to the large scale 
properties of the universe (besides its acceleration) but not vice versa. The standard 
model (and any other model too) provides only some approximation to the reality and 
a slightly modified approximation, which also fits the observations with a satisfactory 
accuracy, would lead to a gravity theory different from general relativity. It has recently 
been found that an attempt to reconstruct the ACDM cosmo logical model as a L = f{R) 
gravity provides the required L = i? — 2A as one of a number of solutions but not as 
the one most probable (i.e. it is not the solution corresponding to attractor points of 
an autonomous system) [29]. In short: from an approximate solution it is impossible to 
reconstruct in a reliable way the correct equation of motion. 
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On the other hand if the exact form of the spacetime metric accurately corresponding 
to the distribution of the cosmic matter were known, the problem of how to reconstruct 
the underlying Lagrangian would be open: the method developed in |28] is specific to 
R-W spacetime and does not work in other cases. We stress that, contrary to a common 
belief, the R-W spacetime is particularly deceptive and unsuitable for recontructing the 
underlying Lagrangian. This spacetime has a high symmetry and is 'flexible' in the 
sense that it contains an arbitrary function, so that it is a solution in any metric gravity 
theory (while Minkowski space is not). Hence for fundamental rather than technical 
reasons it should not be used for the reconstruction. If one believes at all that it is 
possible to recover the Lagrangian from one (empirically found) solution, one should 
apply a solution which does not appear as such in most of gravity theories and is a 
characteristic feature of a possibly narrow class of theories. 

A direct comparison of predictions of a given theory with observations is obstructed 
by the fact that any nonlinear gravity theory may be formulated in infinite number of 
distinct frames and many of them have advantage over the original Jordan frame (in 
which all the theories are initially formulated) in displaying the number of degrees of 
freedom, the particle spectrum and the dynamics of these fields. Since the problem of 
which frame is physical (i.e. consists of directly measurable dynamical variables) still 
remains a matter of a vivid debate and since for this reason any agreement (or dis- 
agreement) of the given theory with the observational data may be criticized, instead 
of attempting to deduce a gravity theory from the data and prior to attempting such 
a confrontation for a chosen theory, one should verify if the theory meets the general 
requirements imposed on a classical field theory. A general L = f{g^u, RafBfiu) metric 
theory has eight degrees of freedom and describes a gravitational triplet consisting of 
the metric, a massive spin-two field and a massive scalar field. For L = f{R) the spin-2 
field disappears. All the L = f{g^u, Rap) theories less differ from each other than it is 
expected in the original Jordan frame since they may be mapped onto general relativity 
(including the two massive fields) and in this sense the latter theory is clearly distin- 
guished as being a universal Hamiltonian image of all these theories. 
How many conditions a gravity theory should satisfy to be regarded a viable one may be 
disputable (e.g. should it have quantization properties better than general relativity?). 
The criterion that a theory have a stable ground state being a maximally symmetric 
spacetime is indisputable. This criterion is investigated in the subsequent Paper II. It 
is shown there that there is an infinity of viable gravity theories. One should therefore 
apply other viability criteria coming from classical field theory. The most natural one 
would be the existence of a proper Newtonian limit. However at present it is difficult to 
formulate this condition in a precise and effective way. And we stress that at the present 
level of knowledge it is impossible to find out a system of selection rules allowing one 
to choose a unique theory out of the whole set of the metric gravity theories. Therefore 
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above all one should provide a deeper physical motivation (different from mere applying 
a number of selection rules and from the wish to account for the cosmic acceleration) 
for choosing a specific Lagrangian rather than any other. In other terms the cosmic 
acceleration should be a new important test for a modified gravity theory but it does 
not provide a way for reconstructing it. 
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